Astronomy & Astrophysics manuscript no. letter 


© ESO 2008 


May 23, 2008 





Missing baryons, bulk flows and the E-mode polarization of the 

Cosmic Microwave Background 

Carlos Hernandez-Monteagudo 1 and Rashid A. Sunyaev 1 ' 2 ' 

1 Max Planck Institut fur Astrophysik, Karl Schwarzschild Str.l, D-85741, Garching bei Miinchen, Germany 
e-mail: chm@mpa-garching . mpg . de 

2 Space Research Institute, Russian Academy of Sciences, Profsoyuznaya 84/32, 1 17997 Moscow, Russia 
e-mail: sunyaev@mpa-garching.mpg . de 

Received ; accepted 

ABSTRACT 

Context. Most of the missing baryons are found in slightly overdense structures like filaments and superclusters, but to date most of 
them have remained hidden to observations. Similarly, the linear cosmological perturbation theory predicts the existence of extended 
bulk flows seeded by gravitational attraction of linear potential wells, but again these remain undetected. 
Aims. Unveil the signature of bulk flows of the missing baryons in upcoming Cosmic Microwave Background (CMB) data. 
Methods. If the peculiar motion of galaxy groups and clusters indeed resembles that of the surrounding baryons, then the kinetic 
Sunyaev-Zel'dovich (kSZ) pattern of those massive halos should be closely correlated to the kSZ pattern of all surrounding electrons. 
Likewise, it should also be correlated to the CMB E-mode polarization field generated via Thomson scattering after reionization. We 
explore the cross-correlation of the kSZ generated in groups and clusters to the all sky E-mode polarization in the context of upcoming 
CMB experiments like Planck, ACT, SPT or APEX. 

Results. We find that this cross-correlation is effectively probing redshifts below z = 3 - 4 (where most of baryons are missing), 
and that it arises in the very large scales (/ < 10). The significance with which this cross-correlation can be measured depends on 
the Poissonian uncertainty associated to the number of halos where the kSZ is measured and on the accuracy of the kSZ estimations 
themselves. Assuming that Planck can provide a cosmic variance limited E-mode polarization map at / < 20 and S/N ~ 1 kSZ 
estimates can be gathered for all clusters more massive than 10 14 M Q , then this cross-correlation should be measured at the 2-3 cr 
level. Further, if an all-sky ACT or SPT type CMB experiment provides kSZ measurements for all halos above 10 13 M o , then the 
cross-correlation signal to noise (S/N) ratio should be at the level of 4-5. A detection of this cross-correlation would provide direct 
and definite evidence of bulk flows and missing baryons simultaneously. 
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, 1 . Introduction that most of the baryons are in a warm phase (WHIM), defined 

• • ■ by the temperature range T e [10 5 , 10 7 ]K. Whether most of 

. ^ I One of the fundamental predictions of the standard cosmologi- this gas is in a diffuse phase or in small, unresolved, collapsed 

V ■ cal model is that baryons (and matter in general) should be mov- objects is sti ll an open question. Recently, |Genova-Santos et al.| 

Sh ' ing in extended bulk flows generated by the gravitational pull of (|2Q05|, |2008|) detected a non-Gaussian cold spot when looking 

^ ! large scale overdense regions. De spite they are potentially useful at the CMB radiation in the direction of Corona Borealis, and 

in terms of cosmological studies jKashlinskv & Atrio-Barandelal discussed its interpretation in the context of the thermal and 

2000j lAtrio-Barandela et al l I2004 |2008l). none of t hose bulk the kinetic Sunyaev-Ze l'dovich effects. The therma l Sunyaev- 

flows have been detected to date dBenson eUd1l2003t) . partially Zel'dovich effect (tSZ, |Sunyaev & Zeldovich| [1980) expresses 

due to the difficulty of detecting baryonic matter in the lo- the distortion that the black body CMB spectrum undergoes due 

cal universe. Indeed, the visible baryonic matter in the var- to Compton scattering on hot elec trons. On the other hand, the 

ious frequency ranges amounts only to one ninth of the to- kinetic Sunyaev-Zel'dovich (kSZ, |Sunyaev & Zeldovich|[l97l 

tal budget of baryons pre dicted by the cosmological model describes the Thomson scattering between CMB photons and 

dFukugita & Peeblesll200l this is the so called missing baryon electrons in which there is no energy exchange (and therefore 

problem). All these baryons must be found in the present it introduces thermal or frequency independent brightness tem- 

Universe, since they have been detected at large (z-l, 100) and perature fluctuations). In the case of Corona Borealis, it was 

intermediate redshift (z ~ 6), by means of Cosmic Microwave suggested that a face on filament could give rise, via either tSZ 

Back ground (CMB) dHinshaw et all l2TjQ8j) and Lyman-ar for- or k SZ, to the pattern seen in the CMB maps . However, it was 

est (iCroft et al.l Il998t iMcDonald et al.1 [2005b observations, shown in|Hernandez-Monteagudo et al.| (|2006|) that the tSZ is not 

However, no conclusive observational evidence has been found a useful tool to track the baryons outside the largest collapsed ha- 

at low redshift. Some indications have been derived from the los : about ~ 80 % of the tSZ luminosity is generated in the most 

obse rvation of ab s orption lines in the direction of quasars massive collapsed structures, and the ~ 50% of baryons located 

(see iNicastro et all d2005albl) and refe rences therein), or a soft in underdense or slightly overdense regions give rise to only ~ 

X-ray excess around clusters, (e.g.. iKull & Bohringeri [l999t 5% of the total tSZ luminosity. The kSZ is most promising, since 

Zap pacosta et al.l2002b . l6erT& Ostriker ( 1999, 2005) concluded it does not require large gas pressure, but relative velocities with 
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Fig. 1. (Left panel, a): Visibility function versus redshift under WMAP5 cosmogony. (Middle panel, b): Thomson optical depth 
versus redshift due to all electrons (solid line) and to those electrons contained in halos more massive than 5 x 1O I2 M (dashed 
line). (Right panel, c, square root of angular power spectra shown, note the p.K units). [Thick solid line:] Total CMB temperature 
anisotropy amplitude. [Thin solid line:] Thomson scattering induced fluctuations during and after reionization. [Thick dashed line:] 
Correlation term of the Thomson scattering induced fluctuations generated in galaxy clusters. [Thin dashed line:] Poisson term 
of the Thomson scattering induced fluctuations in galaxy clusters. [Dot-dashed line:] Cross-correlation between the total E-mode 
polarization (thin solid line) and the correlated part of the kSZ generated in galaxy groups and clusters (thick dashed line). 



respect to the CMB instead (which is itself a prediction of the 
model). However, the drawback is that the kSZ temperature fluc- 
tuations do not depend on the frequency, and therefore they can- 
not be easily distinguished from the intrinsic CMB temperature 
fluctuations. 

In this work we profit the anisotropic nature of the Thomson 
scattering (for which linear polarization is introduced in the 
CMB), and use it as a tool to unveil the bulk flows of the missing 
baryons in the low redshift universe, (z < 2 - 3). For that, we 
need to assume that the peculiar motion of clusters and groups 
are correlated to the velocities of the surrounding matter (just 
as the linear theory predicts), and use measurements on halos 
as a kSZ template for all baryons. In Section (f2| we compute 
the kSZ and the polarization generated during and after reion- 
ization, and the cross-correlation between both. In Section (0 
we compute the kSZ produced by the halo population and the 
signal to noise (S/N) ratio for its cross-correlation to the all- 
sky E-mode polarization. We discuss our results and conclude in 
Section (0]). Throughout this paper, we shall be using a cosmo- 
logical parameter set corresp onding to the WMAP5 + ACDM 
model, dHinshaw et"ai1l2008l) . 



2. The large scale secondary anisotropies after 
reionization 

According to WMAP data, t he first stars must ha ve reionized 
the IGM at around z reio all, dHinshaw et alj|2008l) . Since then, 
CMB photons and electrons have interacted again via Thomson 
scattering. As during recombination, this scattering preserves 
the black body spectrum of the CMB, but however changes the 
direction of the CMB photons, and, in the presence of a non- 
vanishing quadrupole of the CMB intensity, in also induces lin- 
ear polarization. This process is the cause of a partial blurring 
of the anisotropies generated at the surface of last scattering and 
also the generation of new anisotropies at scales comparable to 



the horizon at z — z re io (0 ~ 20° or I ~ 10). (One should also 
consider that, if reionization proceeded first in those most over- 
dense regions hosting the first stars, then it must also have left a 
patchy signature on the CMB small scale anisotropies. However, 
in this work we shall neglect these small angle fluctuations). The 
temperature anisotropies introduced in the CMB intensity were 
caused by the peculiar motions of the scatterers (electrons) with 
respect to the CMB, (i.e., the kSZ effect). A simple expression 
for these temperature anisotropies is given by the following in- 
tegral of the projection of the peculiar velocity along the line of 
sight n: 

— (n) a J dr]A(r])^ J, (1) 

(note that we are neglecting the few percent contribution of 
the polarization terms). As in ISeljak & Zal darriaga dl996l) . in 
this equation r\ is the conformal time (related to the coordinate 
time via dr\ = dt/a(t), with a the cosmological scale factor), 
and T] re io,T]o are its values at reionization and present, respec- 
tively. Dots will denote derivatives with respect to the confor- 
mal time. The (proper) electron peculiar velocity is given by 
\ e (r]), and its Fourier counterpart can be related to the under- 
lying linear density perturbation field by means of the continu- 
ity equation. These peculiar velocities describe the large scale 
bulk flows of matter that are generated by the linear perturba- 
tions in the matter density field. The visibility function is de- 
fined as A(77) = t(t]) exp(-r(?7)), with t(t]) the Thomson op- 
tical depth given by t(i]) = J 10 dr( a(j]')crTn e (ri'). In this ex- 
pression, n e (ri) the background electron number density and cr T 
the Thomson cross section. The Thomson opacity is given by 
t(if) = a(rj)o-jn e (T]). The visibility function A(rf) denotes the 
probability that a CMB photon was last scattered at the epoch 
given by 77. In panel (a) of Figure dTJ it can be seen that the vis- 
ibility function p eaks at recombination redshi ft z - 1, 100 (as 
first predicted by Sunyaev & Zeldovich] dl970h \ although it has 
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a second peak at z — z re io -10, denoting the new Thomson scat- 
tering taking place after the stars reionize the IGM. Note that the 
optical depth behaves differently, increasing steadily upto reion- 
ization redshift (solid line of panel (b) in Figure djjl). 

The anisotropic nature of the Thomson scattering introduces 
linear polarization as long as the CMB intensity has a non- 
vanishing local quadrupole. In CMB studies, it is customary to 
map the standard Q and U polarization St okes parameters into 
the (rotation invariant) E and B modes dZaldarriaga & Seljak 
1 19971: iKam ionkows ki et ai] [l997). which behave with opposite 
parity properties versus spatial reflections. Pure Thomson scat- 
tering should produce no B mode polarization, and hereafter 
this mode will be neglected. Since the evolution of the inten- 
sity and E-mode polarization anisotropics are coupled by the 
sum of the intensity quadrupole and the polarization monopole 
and quadrupole, both quantities will be correlated, and actu- 
ally this cross-correl ation has already been measured by WMAP, 
dHinshaw et al.ll2008l) . 

Using a modifi ed version of a standard Boltz mann code (like, 
e.g., CMBFAST (ISeliak & Zaldarriagal [19961) ) it is straight- 
forward to compute the temperature and E-mode polarization 
anisotropies generated during and after reionization, and the 
cross-correlation betw een them.The result s, computed under 
WMAP5 cosmology dHinshaw et al.l 120081) . are displayed in 
Figure ([TJ;). The total intensity (TT) and E polarization (EE) an- 
gular power spectra are displayed by the thick and thin solid 
lines, respectively, whereas the cross-correlation between the in- 
tensity and E polarization is given by the dotted line (TE). The 
temperature anisotropies induced by the kSZ during and after 
reionization is given by the filled circles, and the generation of 
polarization at the same cosmic epochs is the responsible for the 
low / (/ ~ 5 -8) bump in the EE spectrum. Although the WMAfQ 
mission has barely measured it, the upcoming PlanckQ satellite 
should be able to provide an accurate estimate of the EE angular 
power spectrum (and this shall be our assumption hereafter). 



3. The clusters as probes of bulk flows 

As mentioned in the Introduction, there are two major aspects 
of our understanding of the local Universe that still are pend- 
ing for confirmation: one is the dete ction of the eight ninths o f 
baryons that have not yet been found ( Fukugita & Peebles 2004), 
the another is the measurement of the cosmological bulk flows. 
With the advent of hi gh resolution and high sensitivity CM B ex- 
periments like A CT dFowler & ACT C ollaborationl l2006l) . SPT 
dRuhl et alj2004b or APEX dDobbs et al.ll2006l) . bulk flows have 
prospects to be detected by looking at the kSZ in galaxy clus- 
ters. These structures contain large reservoirs of gas, which 
constitute sources for relatively large values of optical depth, 
(r ~ 10" 4 , 10~ 3 ). Nevertheless, one must have in mind that the 
kSZ effect is not the only physical mechanism to be found in 
clusters. Indeed, given the fact that, in these halos, gas can reach 
temperatures at the level of a few KeV, inverse Compton scatter- 
ing in which hot electrons transfer energy to CMB photons (tSZ) 
becomes of relevance as well. However, the brightness tempera- 
ture fluctuations introduced by this effect cross from negative to 
positive values at a frequency close to 218 GHz. For this reason, 
this frequency becomes an ideal window for searching for the 
kSZ in galaxy clusters. 



1 WMAP URL site: http : //map . gsf c . nasa . gov 

2 Planck's URL site: 

http : / /www . rssd . esa . int/ index . php?pro j ect=planck 



The kSZ effect generated in the gal axy cluster 
population has been studied extensively dPee| [2 006: 
Hernandez-Monteagudo et al. 2006; Bhattacharya & Kosowsky 
2007. 120081 etc). In particular, in iHernandez-Monteagudo et al. 
(120061) an analysis of the correlation properties of the cluster 
peculiar motion was provided. They found that the typical 
correlation length of these peculiar velocities (i.e., the typical 
size of a bulk flow) is typically 20 /z~'Mpc in comoving units. 
This means that a great fraction of the gas that is surrounding 
the clusters is actually comoving with it. This can be envisioned 
as flows of matter (containing clusters, groups and diffuse 
gas) falling into the (future) supercluster's potential well, and 
generating a kSZ dipol e pattern at its cente r as different flows 
merge in that region, (Diaf erio et al.l 12000). This is related to 
the basic assumption of this work: peculiar motion of baryons 
in a given region is well described by the peculiar motion of 
the largest galaxy groups and clusters within the same region. 
When computing the a ll sky kSZ angular power s p ectrum 
generated by clusters, Hernandez-Monteagudo et al. (2006) 
found that although the this could be decomposed in four terms, 
only two of them were giving rise to most of the power: a 
correlation term (generated by the bulk flows) was dominant at 
the very large angular scales, whereas the anisotropy induced 
by the Poisson statistics would dominate at smaller scales. Both 
correlation and Poisson terms are shown by the thick and thin 
dashed lines in Figure dT};). The former has been computed 
using a modified version of the CMBFAST code where the 
Thomson opacity generated in clu sters at a given r\ reads 
dHernandez-Mo nteagudo et al. 2001J: 



r = a(rf)a T \ dM ^zN e (M) b v . 
JM ml „ dM 



dn 



(2) 



The symbol N e {M) - M/(m p fY)Q.b/Q m denotes the number of 
electrons in halo of mass M, with m p the proton mass, fy = 1 - 
Y/2 the Helium fraction correction factor and Q£, Q m the baryon 
and matter density parameters, respecti vely. The halo ma s s func - 
tion dn/dM is approxima ted by that of Sheth & Tormenl (ll999). 
Motivated by iPeell d2006l) . the velocity bias was approximated 
by b v = 1.3 for all halos (although the particular choice of this 
parameter has very little impact on our results, as we shall see 
below). We set M„ mx = 10 16 M , and compared our results for 
different choices of the minimum mass M„„„. In the middle panel 
of Figure ([TJ) the dashed line shows the optical depth generated 
by halos more massive than 5 x 10'2 M : most of the contribu- 
tion (~ 60%) is coming from z < 3, and the contribution from 
(the more massive) galaxy groups and clusters is shifted even 
further into the low z range. On the other hand, the Poisson term 
can easily be found to be 



C 



kSZ£ 



p(,„„„) rM„ ulx dn 

dnn 2 dM 

~L, Jm, u ,„ dM 



o- T N e {M)b v o- v {M,rj) 



Ang 



(V) 



x [1 + ^,1. 



(3) 



where aUn#(?7) denotes angular diameter distance and cr v (M, rj) 
is the time dependent rms line of sight peculiar velocity of an 
object of mass M, 



o 2 v {M,ri) = \\H( U ) 



dz 



)7 



(4) 
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Fig. 2. The three panels display the total S/N ratio for the correlation of the E-mode polarization with the low redshift (z < 3 -4) kSZ 
generated by bulk flows in galaxy groups and clusters. Horizontal axes display the decimal logarithm of the individual group/cluster 
kSZ measurement S/N ratio (£), whereas vertical axes display the decimal logarithm of the minimum mass of groups/clusters 
considered. Due to Poisson/shot noise, a low choice for Zmin makes the high total S/N regions shrink, and this also happens if z mi „ 
is chosen too big and too much information is unused. The optimal minimum redshift is around z„„„ ^ 0.05, and the maximum S/N 
ratio achievable for perfect E and kSZ surveys is ~ 5.8. Note that we are assuming /,a v = 1. 



The present linear matter power spectrum is given by P m (k), 
the linear density growth factor by Ds, the Hubble parame- 
ter by H{rf), and the Fourier window function of a top hat fil- 
ter corresponding to an objec t of m ass M by W(kR(M)), (see 
iHerna ndez-Monteagudo ^tail d2006l) for details). Note that in 
equation (0, £ denotes the signal to noise ratio of the kSZ de- 
tection in a single cluster/group. As noise we regard any signal 
present at the cluster position that is not due to the correlated 
component of the kSZ (e.g., kSZ caused by thermal velocities, 
tSZ residuals, intrinsic CMB, radio source emission, IR source 
emission, instrumental noise, etc). Note as well the low redshift 
limit given by z m! „ in equation (|3). Therefore, this equation ac- 
counts for both the Poissonian statistics and the error in the kSZ 
recovery at cluster positions. Only in the case £ — > oo we recover 
the pure Poissonian term. 

We are interested in the correlated component of the clus- 
ter kSZ, i.e., in the velocity component that should be comoving 
with the surrounding matter. If this is indeed the case, then the 
pattern of the kSZ in clusters must be very similar to that of the 
kSZ generated by all electrons (filled circles in Figure ([TJ), and 
therefore it must be closely correlated to the E mode of the polar- 
ization anisotropies generated during and after reionization (i.e., 
the low I bump). The actual amplitude of this cross-correlation 
will depend on the redshift range that a given set of halos is prob- 
ing, and the signal to noise (S/N) ratio of this cross-correlation 
for a particular multipole I is provided by 

/£_ \ 2 _ & + I) f sky (cf M f 5 

\N II qEE ^jkSZ.cl + gkSZ,^ + lgkSZ,cl-E^ 2 ' 

The symbol f s k y refers to the fraction of the sky covered by an 
experiment, and c* SZc/ and cf z ^ denote, respectively, the an- 
gular power spectrum of the correlated part and the noise of the 
kSZ generated in galaxy clusters. The cross-correlation between 
the kSZ in halos and the E-mode is given by C* sz ' e ' _£ . If / s t v = 1, 
the total S/N is given by the sum of the squares of these quan- 



r Til/2 

tities, i.e., S/N = 1 2/(5 /N)f I . It turns out that the actual am- 
plitude of this S/N ratio is strongly dependent on the relative 
amplitude of the Poisson term with respect to the correlation 
term (cancelling any dependence on parameters like the veloc- 
ity bias £>,,). Due to the rj 2 /d^ °c l/rf dependence in equation 
©, the Poisson term is less important for higher values of the 
minimum redshift considered Zmin- Further, it depends linearly 
with the number density of objects, as opposed to the correlation 
term, which depends on the number density squared. Therefore, 
for smaller values of the threshold mass M m ,„ the relative weight 
of the Poisson term versus the correlation one should diminish. 
This is shown in Figure d2j. It is worth to remark as well that 
since this correlation arises around halos at low and moderate 
redshifts, the matter behind this signal is different from that gen- 
erating the and cross-correlation between CMB anisotropies and 
HI 21cm observations. 



4. Discussion and Conclusions 

The horizontal axis of Figure (fJJ corresponds to the S/N ratio 
to which the kSZ can be determined in a single halo, i.e., £ in 
equation d3). The number of potential contaminants when mea- 
suring the kSZ in a halo is relatively high: the tSZ has a definite 
spectral dependence, and should be close to zero at v = 218 
GHz, but however asymmetries in the freque ncy response of 
the ex periment together with relativistic effects (lltoh & Nozawal 
2004) will inevitably leave some resid uals. However, as shown 
in lHernandez-M onteagu doet al.l d2006l) . its relative weight to the 
kSZ decreases for the most numerous low mass halos, for which 
the relativistic tSZ corrections are also less imp ortant. Internal 
motion s within the halos (which, as found by iDiaferio et al.l 
(2005) are only relevant for slow halos) should average out 
when integrating over the total halo solid angle. The uncorre- 
lated component of the kSZ could be significant in those sys- 
tems where the velocity of halos are thermalized; however this 
almost never happens for the galaxy groups and clusters being 
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targeted here, since supercl usters are structures that are not yet 
relaxed. Nevertheless, iPeell ([2006 ) found that, in his cosmolog- 
ical simulations, the peculiar velocities of haloes showed cor- 
relation properties that differed substantially from linear theory 
predictions, even at high redshift. If this is indeed the case, then 
the issue would be whether the baryons surrounding halos would 
have their peculiar velocities affected in the same way or not. We 
shall address this issue in future work. The impact of IR and ra- 
dio sources is a priori significant, provided that these sources 
are correlated to the halos where the kSZ is produced. However, 
multifrequency observations should provide a good handle on 
their spectral indexes, and therefore their contribution at the 218 
GHz channel should be characterized down to an accuracy level 
to be realistically estimated only when real data becomes avail- 
able. The intrinsic CMB generated at the last scattering sur- 
face is anoth er potential source of contamination. As shown in 
Figure (6) of Hern andez-Monteagudo et ail d2006l) . its residuals 
remain in the few pK level if the cluster/group size is not larger 
than 2-3 arcmins, which should be the case for most of clusters 
and groups in the sky. Therefore this source of contamination is 
likely to be less relevant than the rest. 

The importance of each of these contaminants will determine 
which range of the £ horizontal axis in Figure (f2]i one should look 
at. Due to the largest amplitude of Poissonian fluctuations at low 
values of z mlH , it is convenient not to consider too nearby clusters 
and groups. From this figure it is possible to find at which values 
of £ the S/N ratio starts being dominated by Cosmic Variance 
(horizontal plateau of the contour levels at high £). Likewise, 
the steep vertical descent of these contours at low £ reflects the 
limit where S/N is exclusively limited by the accuracy of the 
kSZ measurements. The requirements for a convincing detection 
of the signature of the local baryon bulk flows (say S /N > 3) in- 
volve either a large value of £ or a low value of M„„„, but never- 
theless these are not unrealistic. For instance, the Planck satellite 
should provide a cleaned CMB map in a very large fraction of the 
sky, and with it, estimations of the kSZ at the cluster and group 
positions. These estimates should mostly be limited by the beam 
size and the un-substracted point source emission. By looking 
at the Compton distortion y map (obtained from multifrequency 
observations of the tSZ), one should be able to find the position 
of not only the largest clusters, but also of those objects that lie 
closer to the confusion limit, and whose signal can be picked up 
statistically via the cross-correlation. If this y-map provides the 
position and the cleaned CMB map the kSZ estimates (of at least 
£ ~ 1) for all clusters more massive than 1O 14 M , then Planck 
alone would already yield a total S/N between 2 and 3. Likewise, 
upcoming Dark Energy surveys like Pan- STARRS0, DES n or 
PAU-BAC0, or X-ray surveys like eROSITA0 should probe the 
relevant redshift range, and provide independently catalogs of 
group and cluster candidates. At the same time, high resolution 
CMB experiments like ACT, SPT or APEX would contribute 
with more accurate kSZ estimates. Furthermore, future B-mode 
polarization experiments should cover large fractions of the sky 
with very high sensitivity, and provide better understanding of 
the E-mode of the CMB polarization. One must remark that the 
maximum achievable value of the total S/N is relatively small 
(around 5.8), and therefore a large sky coverage is required, 
since most of the signal is arising at Z < 10. Along these lines, 
we must stress again that a cosmic variance limited EE power 



spectrum has been assumed, at least in those same large angu- 
lar scales. Very soon, the Planck mission should be able to tell 
whether this is a sensible assumption or not. 

Acknowledgements. C.H.M. acknowledges useful conversations with 
D.Spergel. 
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